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The eigenenergies e'^' (m; {n^, /i}) of a system of A'' identical particles with a mass m are func- 
tions of the various radial quantum numbers Ui and orbital quantum numbers k. Approximations 
E^^\m;Q) of these eigenenergies, depending on a principal quantum number Q{{ni,li}), can be 
obtained in the framework of the auxiliary field method. We demonstrate the existence of numer- 
ous exact duality relations linking quantities E^^\m;Q) and E^^\m'; Q') for various forms of the 
potentials (independent of m and A^) and for both nonrelativistic and semirelativistic kinematics. 
As the approximations computed with the auxiliary field method can be very close to the exact 
results, we show with several examples that these duality relations still hold, with sometimes a good 
accuracy, for the exact eigenenergies e^'^\m; {ni,li}). 
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; I. INTRODUCTION 

. The auxiliary field method (AFM) is a very poviferful method to obtain analytical expressions for the eigenvalues 
of one, two and many-body systems. At the beginning, it was introduced to deal, in a simple way, with the worrying 
square root operator appearing in many relativistic theories [H-Q- But we showed, in a series of papers (5l-[llj|. that it 
^ ' can be applied as well with great success in many physical situations. In this work, we focus mainly on the properties 
^ I of the eigenenergies of iV-body systems and the relations linking them. So, we will recall only the basic ingredients 
^H. necessary for the understanding of the subject treated here, and we refer the reader to our review paper [ 12{ for an 
exhaustive overview of the method and its applications. Many detailed discussions on the AFM properties for A^-body 
systems can also be found in 

A great advantage of the AFM lies in its great simplicity and its ability to treat on equal footing the ground state 
and the various excited states. An analytical expression for these eigenenergies (the mass M for a semirelativistic 
^vq ' approach or the binding energies E for a nonrelativistic one) is particularly interesting from two points of view: 



m 
O 



• It allows a very fast calculation in terms of the parameters and thus is very well suited if such calculations are a 
part of a chi-square determination of these parameters, a procedure which needs a lot of calls to the calculation 
of the energies. 

• It provides the behavior of the eigenenergies as functions of the various parameters and quantum numbers. 



This last item can be invoked to search for relationships between the eigenenergies of different systems with different 
, y , sets of parameters: we call these relationships "duality relations" . They form the main subject of this paper. 
rN In essence, the AFM replaces a problem which is not solvable, for example because of a complicated potential V{r), 

by another one which can be treated analytically, for example by the use of a more convenient potential P{r). In so 
" ■ " doing, it is necessary to introduce an auxiliary field v which is an operator. The original Hamiltonian H is replaced 
by a new simpler Hamiltonian H{v)^ called the AFM Hamiltonian. If this auxiliary field is chosen as in order to 
extremize the AFM Hamiltonian, it can be shown that the AFM Hamiltonian coincides with the original Hamiltonian: 
i7(j)o) — H . Thus, both formulations are completely equivalent. 

The approximation lies in the fact that the auxiliary field is considered no longer as an operator, but as a real 
constant v. The value vq which extremizes the eigenenergy E{v) of the AFM Hamiltonian is then inserted in the 
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corresponding expression to give an approximate value E{vq) of the exact eigenenergy E. An approximate state for 
the corresponding eigenstate can also be obtained. The quality of this approximation has been studied and discussed 
in detail in the papers mentioned above. 

In this paper, we are concerned with systems composed of N particles interacting via one-body potentials Ui and 
two-body potentials Vij, and moving with a nonrelativistic or a semirelativistic kinetic energy. Although the AFM 
can be invoked to treat this problem in such a general form, it is manageable in practice only if the particles are 
identical. This implies that they all have the same mass m, that the form of the one-body potentials is the same for 
all particles Ui = U, and that the form of the two-body potentials is the same for all pairs of particles Vij = V. In 
consequence, the Hamiltonian under consideration in this work has the following form 

N N N 

H = Y.^pUm^ + Y,U{\n-R\)+ ^ nh-r,|). (1) 

i—l i—1 i<j—l 

In this expression is the position operator for particle i, pi its conjugate momentum and R the position of the 
center of mass of the N particles. 

The only case which is entirely solvable analytically is a nonrelativistic system with quadratic potentials U{r) = 
V{r) = r^. With the natural choice P{r) = r^, the AFM needs two auxiliary functions (the prime symbol stands for 
the derivative with respect to the argument) K{r) = U'{r)/P'{r) ^ U'{r)/{2r) and L{r) = V'{r)/P'{r) = V'{r)/{2r), 
and the introduction of three auxiliary fields: /i to replace the relativistic kinetic energy by a nonrelativistic one, v to 
replace the potential U{r) by the potential P{r) ~ r"^ and p to replace the potential V{r) by the potential P{r) — . 
The minimization procedure giving the best values /io, '^Oi Po leads to a system of 3 coupled non linear equations (see 
[lllj). 

Introducing the new variable 



Xo = ^/2iI^^(^J^^+N]^, (2) 

it appears that the previous system can be recast under the form of a single transcendental equation which looks 
like 



NXo / \ V - 1)^0 



(3) 



while the corresponding AFM eigenmass is given by 



where the number of pairs 



C. ^ M (5) 



has been introduced for convenience. 

In these expressions, Q denotes the principal quantum number. Concerning this point, some comments are in 
order. The system being described by {N — 1) Jacobi internal coordinates, the most general excited state depends 
on {N — 1) radial quantum numbers rii and {N — 1) orbital quantum numbers li, as well as intermediate coupling 
quantum numbers which are not considered here. The AFM relying on a quadratic potential, the principal quantum 
number resulting from the AFM treatment is naturally 

N-l 

^(2n, + /,) + -(iV-l). (6) 



In particular, the ground state for a A''-boson system is just Q — 3{N ~ l)/2, while the ground state of a A^-fermion 
system is much more involved and needs the introduction of the Fermi level 

In a two-body system, Q reduces to 2rt -I- ^ -I- 3/2. But it is possible in this case to use other forms for P{r) leading 
to other expressions for Q. For instance, Q = n + I + 1 with P{r) = — 1/r @ and Q = 2(— q;„/3)'^/^ for S-states with 
P(r) = r, where a„ is the (n -I- 1)**^ zero of the Airy function Ai fT^]. Moreover, it has been shown (for A^ = 2 in 



Q and for arbitrary A'' in that a much better approximation of the exact energies can be obtained with a shght 
modification of the principal quantum number. A particularly simple form which works quite well is given by 



In the following, we consider the principal quantum number Q as a variable by its own in all the computations, and 
only at the very end of the calculation we replace this undetermined variable by the most precise one Q = Q({7ii, li}) 
as in ^ or in ([T]). 

It is remarkable that, for such a complicated and sophisticated system, the only difhculty to get the AFM approx- 
imate solution needs to solve the transcendental equation ([3]), a quite easy procedure numerically. Moreover, in a 
number of interesting cases, this solution is analytical so that the mass itself (given by (0])) is also analytical. 

In general, the potential U{r) (or V{r)) depends on physical parameters {n, T2, . . . , Tp} = {r} and should be noted 
more precisely J7({r}; r). The other parameters of the problem are the number of particles N and the mass m of the 
particles. Lastly, we want to have a description of the whole spectrum, so that the principal quantum number Q also 
enters the game. In consequence, a complete notation for the eigenmasses would be Af'^^^({T}; m, Q). 

Let us stress now a very important point: for realistic potentials, the parameters {r} could depend on N or/and 
m. We do not consider such a particular behavior here. Thus, in this paper, we assume that the parameters of the 
potentials are independent of N and m. This means that the m dependence of H is only through the kinetic 
energy and its N dependence through the numbers of terms in the summation. In this paper, we suppose that the 
potentials are given once for all and consider that their form and parameters do not vary for all studied systems. 
Therefore we are finally interested only by the N, m, Q dependence of the eigenmasses and use the simplified notation 



The two-body problem, N = 2, deserves a special treatment. In this case R ~ {ri+r2)/2 so that |ri — ii| — \r2 — R\ 
= (l/2)|ri — r2\. Thus, U{\ri — R\) + U{\ri — R\) — 2f/(r/2) where r is the inter-distance between the particles. 
We can then write 2[/(r/2) = V{r). In consequence, the distinction between one-body and two-body potentials 
becomes irrelevant and it is sufficient to deal with only one type of potential and, depending on the situation under 
consideration, we will choose the most convenient one for our purpose. 

Moreover, the semirelativistic kinetic energy writes (p is the conjugate momentum of the relative distance r) 
T — a^yp^ + rnP , with a = 1 for one-body systems and a — 2 for two-body systems. However, we found very 
convenient to consider the a parameter as a free real parameter, because the eigenenergy can be calculated as well as 
a function of cr. This possibility leads to interesting properties. Thus, eigenmasses of the systems with the previous 
form of the kinetic energy will be noted M'^^\a,m,Q) while the notation AI^^\m,Q) is reserved to the natural 
definition M(^=^^m,Q) = M(^)(2,m,Q). 

What we call a "duality relation" is just a relation between M'-^\m,Q), M'^P\m' ,Q') and M(^)((t, m", Q"). The 
search for such duality relations is the subject of this paper. Of course, these duality relations are exact only for the 
AFM approximation of a given problem. However, since this approximation may be, under many circumstances, a 
rather accurate result of the exact eigenvalues, our hope is that these duality relations remain approximately valid for 
the exact results and thus represent a rough view of reality. In particular, starting with the expression for a two-body 
system, whose calculation is rather easy, one can hope that the main trends for a A'^-body system are reachable as the 
consequence of the duality relations. 

In the next section, we discuss duality relations in the general case, while in the 3rd and 4th sections, we investigate 
the particular cases of ultrarelativistic systems and nonrelativistic systems, respectively. Their interconnexion is 
investigated in detail in the 5th section. In the 6th section, applications are given while conclusions are drawn in 
the last section. Some checks of the duality relations obtained are performed with explicit AFM solutions in the 
Appendix. 



N-l 




(7) 



i=l 



MW(m,g). 



II. GENERAL CASE 



A. Compact form 



In this section, we show that it is possible to solve formally the transcendental equation ([3]) and express the AFM 
mass (U) under a more compact form. In order to do that, we switch from the Xq variable to the y variable defined 
by 
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With this definition, it is easy to calculate the following quantities: 
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TV 2/^ ^' V miV^' Y(iV-l)Xo 777VC;7' 

Let us introduce now, instead of potentials U and V , the related potentials W ^ Y and Z by 



777 



N 



V ( — ^=a; ) = NY{x), (10) 



777VC7V 

It is easy to prove that the transcendental equation ^ can be recast under the reduced form 



y^ = y'Vl+¥z'iy). (11) 
The new step of our procedure is the introduction of an inverse function A, which is defined through the relation 



A{x^y/l + x^Z'{x))^x. (12) 
One can thus express formally the value of y in term of the A function as 

The eigenmass follows from (|4]) introducing rather the y variable. Explicitly, 



M = Nm^^^—^^ + NCn Z(y). (14) 

y 



The last step of our treatment is the introduction of the B function as 



Bix) = ^^^l + Aixy+ZiAix)). (15) 

One sees that the eigenmass is expressed simply in a very compact form 

M^^\m,Q)^NCNB(^^y (16) 

It is astonishing that the eigenmass of such a complicated system, even in its AFM approximation, can be expressed 
in such a compact form in term of a single function B{x). However, this enthusiastic conclusion must be moderated 
by the following remarks: 

• The B function is expressed through the A function which, itself, needs the inversion of a function which 
contains a square root and a potential with a complicated form; this can be a tremendous task. The effort of 
solving the transcendental equation ([3]) is replaced by the inversion procedure (|12p . which is formally as much 
as complicated. 



• As seen from (|10p . the A function, and hence the B function, depends not only on the considered potentials - 
an unavoidable ingredient - but also on N, m and Q, and should be noted more precisely A{N ,m^Q\x). This 
means that this function depends not only on the considered systems (through the potentials, and N and 777 
variables) but also on the particular excitation state Q. 

In consequence of these arguments, it appears that this formulation is of no practical use and do not lead to duality 
relations (at least we did not succeed to find one!). However, we gave this demonstration because, independently of 
the formal esthetical presentation, it represents the prototype of many similar approaches that will be extensively 
used in the following sections. 
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B. One-body interaction only 



In absence of two-body interaction, V = 0, simplifications occur and one can draw interesting conclusions. Instead 
of the y variable, let us use rather the tq variable (which is called the mean radius) defined by 



Q 

mN 



y 



NXo 



(17) 



With this new variable, the transcendental equation (jlip and the mass expression are replaced by the following 
ones 



Q 

N 



= rWl + N^mro/Q)^ C/'(ro), 



Af(^)(TO,g) ^ N 



Q 



^l + N^{mro/Q)^ + C/(ro) 



(18) 
(19) 



In this case, we do not want to get a compact closed formula (one thus has to derive tq from the first equation giving 
ro(iV, m, Q) and inserting it in the second equation to get M'^^\m^Q)). We prefer to express in the same form the 
case of a system of p particles of mass m! interacting with the same one-body potential U (r) (don't forget that it 
is N and m independent) in a state of excitation Q' . We thus have 



Q' 



M'^P\m',Q') = p 



Now, let us choose the parameters m' , Q' as 



Q' 



N 



Q. 



(20) 
(21) 

(22) 



We suppose that the function 
a^x^U^x) is monotonic in x for physical values of x and a (here a = mN/Q); this is always the 



Inserting these values in ([^0]) gives exactly equation ([T51) with sq in place of rg 
f{a;x) = x^VT 

case in practice. Since we have /(a;so) = /(a;ro), this implies sq = vq. Inserting these values of m', Q' and sq in 
(PTjl . it appears that M^p'' = {p/N)M^^\ We are thus led to the following duality relation between the A^-body and 
p-body systems governed by the same one-body potential 



p \ N / 



(23) 



It means that the spectrum of a iV-body system is the same as the spectrum of a p-body system (with the same particle 
mass and the same form of the one-body potential) provided we look at different excitation states. Formula tells 
that the energy is equally divided on all the particles. This is formally exact for arbitrary values of the principal 
quantum number Q; however if one relies on the value Q{{ni,li}) as in (jS]), there is no certainty that pQ/N appears 
still under the form (jG]). Nevertheless, as we already stressed, we must make the calculation of the p-body system in 
term of an informal Q and not trying to fit the spectrum of the system, then make the replacement Q — > NQ /p to get 
formally M^^\m,Q) and, lastly, use in this expression the value of Q = Qilniji}) that fits at best the spectrum. 
These comments are important; in the following we must keep them in mind for all duality relations. 

Let us consider a 2-body system, which is the most simple system to be studied. To apply the previous formula, 
one must take some caution with the definition of the potential since V{r) = 2U{r/2). Thus, if M<-'^\m,Q) is the 
eigenenergy resulting from the Hamiltonian 2y^p'^ + rn? 4- ^(r) (which has the form of a spinless Salpeter equation 
with a two-body potential V), the calculation for M^^^ (to, Q) must be done with one-body potential U {x) — V{2x) /2. 
This means that the determination of the mass of the TV-body system results rather from the equations 



MW(TO,g) 



N 



4g 



Vl + ^2(m.so/(2g))2 + T/(so) 



(24) 
(25) 
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It could be interesting to deal with the quantity M*^*^ (cr, m' , Q'), eigenvalue of the Hamiltonian a ^/p^ + rn/^ + V (r) . 
In this case, M^^^ is given by the following set of equations 

aQ' - rl^l + im'ro/QT V'{ro), (26) 
M(^)(a,m',Q') = ^ + K^o/Q')' + ^M- (27) 

Putting the value m' — 2m /a and Q' — AQ/{aN) into the previous equations, one sees that we recover ([M)) provided 
that To = So and, then, Af*^"^ = {2/N)M^^\ Therefore, one has the duality relation 

MW(m, Q) = |m(^) (^a, ^m, . (28) 

This relation is particularly interesting because it is valid whatever the value of cr. In particular, the special value 
a = 2 gives the hnk between M^^) and M^^), as in ((23)) . 
Choosing the value a — 4/N leads to the duality relation 

MW(m,Q) = ^Af(^) (^l,^m,g) . (29) 

In this case, there is a one-to-one correspondence of the spectrum at the price of considering systems with different 
particle masses. The conclusion is that, if A/^*^ or M^^^ can be evaluated analytically, the same property holds for 



C. Two-body interaction only 



We consider now the case of systems governed only by two-body forces, U — 0. One can model our reasoning on 
the demonstrations given in the previous section. In this case the natural variable is the mean radius 



'"0 



Q 



y 



2Q 



{N-1)Xq- 



(30) 



With this new variable, the transcendental equation (|lip and the mass expression are replaced by the following 
ones 



2Q 



M^'^\m,Q) = Cn 



rl^l + CN{mr^/QY ^ (ro) 
2Q 1 



{mrQ/QY + y(ro) 



(31) 
(32) 



Exactly in the same way as the proof given before, one arrives at the following duality relation between the iV-body 
and the p-body systems with particles interacting via the same two-body interaction 



p — 1 P ~ 1 
-m. 



(33) 



In this case, the spectrum of the iV-body system is the same as the p-body system (with the same two-body potential) 
provided we consider different particle masses in both situations and for other excitation states. 
The hnk between Af and A4"(*) now looks hke 



Q 



(34) 



Here again, this relation is valid whatever the value of cr. In particular, the special value cr = 2 gives the link between 
Af(^) and Af^^^, as resuhing from with p = 2. 

Alternatively, one can use for instance this freedom on cr to choose the same mass in both systems (this was not 
possible with the expression A^^^^); indeed, with the choice a = 2/{N — 1), one has 



A^W(m,Q) = C7jvAf(") 



(35) 
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Choosing a = '2,/{{N — \)\/Cm) allows to express the duality relation keeping the same principal number in both 
systems 



M^^\m,Q) = CnM^ 



{N-1)VC^ 



(36) 



D. Link between one and two-body interaction 



We consider the case of a A^-body system with particles of mass m and interacting with a two-body interaction 
V{x) only. The corresponding mass for a state of excitation number Q, M^^\m, Q), is given by the set of equations 
(1511) and ([5^ . Now we consider the case of a A^-body system with particles of mass m' and interacting with one- 
body interaction of the form U{x) = cV{x) only (c being a constant). The comparison for the mass of both 
systems may present some interest: for example in 3-quark systems, the case U{x) — ax corresponds to the so-called 
Y confinement potential with the junction point placed at the center of mass, while the case V(x) = ax/2 (in this 
case c — 1/2) corresponds to the alternative description of confinement named the A potential [13j. In this case, the 
corresponding mass for a state of excitation number Q', M^^''{m', Q') is given by the set of equations (|20| and ((2T|) 
with p = N. 

If one chooses the mass and principal quantum numbers fulfilling the following relations 



2cQ 



91. 

TV' 



(37) 



it is easy to show that ([20)1 is transformed into (I^Tj) while M^^^/{cN) is transformed into M^^'' /Cn- Therefore we 
arrive at the following interesting duality relation 



2c 



2c Ac^/C^ 
iV^""' (iV- 1)2 



(38) 



III. ULTRARELATIVISTIC LIMIT 

The case of ultrarelativistic systems, characterized by a vanishing mass m — 0, presents some very specific and 
interesting features. For instance, the case of systems composed of gluons or/and light quarks can be well represented 
in this scheme [l3|- The Hamiltonian of the system is then 



N N N 

H = J2^^ + J2u{\r,-R\)+ Yl V{\r,-r,\). 



(39) 



Indeed, the formulation is simpler for this particular situation. Putting the value to = in ([3]) and (Hj), one gets a 
new set of equations. The transcendental equation looks like 




K 



Q 



while the corresponding AFM eigenmass is given by 



NL 



2Q 



iN-l)Xo 



(40) 



(41) 



The eigenmass, depending now only on N and Q, will be noted Mu^\q) (the index u stands for "ultrarelativistic" 
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A. Presence of one and two-body potentials 

One can mimic there the approach foUowed in section Hi Al Let us introduce, instead of Xq, the new variable y 

y = M (42) 

V Aq 

where a is any real parameter. 

Let us define now, instead of potentials U and V, the related potentials W, Y and Z by 

u(^=x] = CnW{x), 



/aN 



Z{x) = W{x) + Y{x). 

It is easy to prove that the transcendental equation (|40p can be recast under the reduced form: 



= y'Z\y). (44) 



NCn 

As in Sect. Ill A[ we introduce A„, the inverse function of y'^Z'{y), 

Au{x^Z'{x)) ^ X. (45) 
One can thus express formally the value of y in term of the function as 

The eigenmass follows from (|4ip introducing rather the y variable. Explicitly 

^ Q 



Lastly, one defines the function i3„ as 



NCNy 



Z{y) 



(47) 



B^{x) = -^ + Z{A^{x)). (48) 



The eigenmass (|T7)) is expressed very simply as 



MW(Q)=iVC^S„(^^^j. (49) 

As in Sect. Ill Al we are able to express the ultrarelativistic eigenmass in a very compact form through a single 
function i3„. Nevertheless, in contrast to Sect. Ill Al we have in this case more sympathetic features: 

• The inversion procedure (1451) is much easier than the corresponding one (|12p , because of the disappearance of 
the square root. 

• The parameter a is free, so that we can choose it in the most convenient way. It appears in the argument of the 
Bu function, but this function itself, as can be seen from the definitions (021), depends on a, so that the final 
result Af„ is a-independent. 

• The potentials W, Y, Z, and, thus, the function Bu does depend on U, V, a and N, but not on Q. Therefore 
the only dependence in Q of the final result is through the argument of the i?„ function. 

As a consequence of the last item, for a given system {U, V and N are given), if there is some degeneracy due to 
Q (for example using the harmonic oscillator expression ([6])), this degeneracy persists in the AFM expression (|49)) of 

the spectrum. Moreover, if it is possible to obtain an analytical expression for Bu{x), the behavior of m2^\q) as 
function of Q directly follows from (HH)) . 
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B. One-body interaction only 



In absence of two-body interaction, V — 0, further simplifications occur and lead to very interesting conclusions. Let 
us choose the special value a = 1/N in the previous approach. One has ^/oN = 1 and Z{x) = U{x)/Cn- Introducing 
the mean radius tq = y, the transcendental equation (|44[) reduces to 

^=rlU'{ro). (50) 

Let us denote C{x) the inverse function of x'^U'{x), 

C{x^U'{x)) = X. (51) 



Therefore, the mean radius is given by 



Lastly, introducing the function F{x) by 



the mass of the system is given by 



F{x) = -^^+U{C{x)), 



(52) 
(53) 
(54) 



The great advantage of this formulation with respect to the approach of Sect. IIII Al is that the C and F functions 
are universal in the sense that it depends on the form of the U potential, but it is independent of the system and of 
the excitation number (independent of N, Q). For two different potentials. The C and F functions, although denoted 
by the same label, are indeed different. 

From expression Ij54p . one deduces immediately the duality relation 

mW(Q) = ^A/(^)(^q), (55) 

which is the special case of ((23|) with m = 0. 

Let us consider the Hamiltonian H = a^/j^ + V{r), whose eigenmass is noted Mt\a,Q). It is easy to show that 

Mi'\a,Q)^F{aQ), (56) 

with the function F as in ([55]) and C{x) as in ([5T|) but calculated with the potential V{x). 

The universality of F{x) does not bring fundamental new results. In practice, all the properties derived in Sect. Ill SI 
can be invoked with the special choice to = to get simplified formulations. For example, duality relation (1^51) writes 
now 

Mr^(Q) = fM(^^(-,^Q), (57) 

valid whatever the value of a. Putting ct = 2 in this relation, one recovers (j55p with p = 2. Using a — 4/A^, one 
obtains the following interesting duality relation 

MW(Q) = fMl^)(|,Q), (58) 
making a one-to-one correspondence in the spectrum of both systems. 
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C. Two-body interaction only 

In absence of one-body interaction, U — 0, we have also interesting conclusions. Let us choose the special value 
a — 2/{N— 1) in the previous approach. One has V aN ~ N/\/Cn and Z{x) = V{x)/N. Introducing the mean radius 
To y, the transcendental equation (|44)) reduces to 

rlV'ir,). (59) 



{N-l)^ 

Defining the C{x) and F{x) function as in ((5T|) and (|53)) calculated with the potential V{x)^ one has 

-^Kra^)' '''' 

Mi^\Q) ^ C^F ( . (61) 



As the previous case studied in Sect. IIII Bl the C and F functions are universal. Let us note that (|5T|) with N — 2 



coincides with (l56l) for a — 2, as expected. 

From expression (j6ip . one deduces immediately the duality relation 



which is the special case of ([55]) with m = 0. 

Application of (j34l) with to = provides the duality relation 

«<"'«) = C„Mi..(.,-^^^^-i^0), (63) 

valid whatever the value of a. Putting tr = 2 in this relation, one recovers (l62|) with p — 2. Using a = 2/ {{N —1)\/Cn), 
one obtains the following interesting duality relation 



M^HQ) = Mli'^^ ( Q ] , (64) 



making a one-to-one correspondence in the spectrum of both systems. 



D. Link between one and two-body interaction 



With the same conditions of application studied in section Hi Dl the interesting relation ([55)1 looks even simpler 

2c " V(^-l) 



^eHQ) = ^A?r^ ( ] ■ (65) 



IV. NONRELATIVISTIC LIMIT 



Another interesting limit of the theory is the nonrelativistic one, valid when the mass of the particles is large 
compared to the mean potential. In this case, the considered Hamiltonian is simply 

^ = E£ + Ef^(i^'-^i)+ E n\n^r,\). (66) 

z— 1 i—1 i<ij—l 

Instead of dealing with the total mass it is better to consider the binding energy obtained by removing the 

total rest mass: approximation E^^\m,Q) of the binding energy is given by the 

following equations. 



= 2m 
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and 

We will see that, in this nonrelativistic limit, one obtains additional interesting properties. 

A. Presence of one and two-body potentials 

Let us introduce the same parameter y as in and the same Z{x) function as in (|43p . The transcendental 
equation looks like 

y^'Z'iy). (69) 



tuNCn 

We introduce An, the inverse function of y^Z'{y) 

A„{x^Z'{x)) =x. (70) 
The value of y in term of the function An is expressed as 

V-aJ^). (71) 



(72) 



ttiNCn 

The eigenenergy follows from (|68p introducing rather the y variable. Explicitly 

aQ2 



E^^\m,Q) ^ NCn 
Lastly, one defines the function _B„ as 



2mNCNy^ 



z{y) 



B^{x) = ,^^-^ + Z{A^{x)). (73) 
The binding energy (|72p is expressed very simply as 

i?(^)(m,Q) = NC^Bn (^$A ■ (74) 

As in the case of ultrarelativistic limit, the function i?„ does depend on U , V, a and N, but not on m and Q. 
Thus the conclusion concerning the degeneracy of the spectrum remains valid. Moreover, if one knows an analytical 
expression of Bn{x) then (l74l) provides the behavior of E^^^ {m,Q) as a function of m and Q. 

In the case of a nonrelativistic system, an additional property appears. Remarking that iJ^^^ depends only on the 
ratio Q^/m, one has the general duality relation 

£;W(to,q) = £;W(/32,„,/3g), (75) 

valid for any value of the real parameter /3. This relation is very strong because it is completely general, even if the 
system is governed by the presence of both one and two-body interactions. 

B. One-body interaction only 

As in Sect. IIIIBl the absence of two-body interaction, V — 0, implies further simplifications. Let us choose the 
special value a = 1/N in the previous approach. One has V aN — 1 and Z(x) = U{x)/Cn- Introducing the mean 
radius tq — y, the transcendental equation (|69p reduces to 



(76) 
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Let us denote -0(2;) the inverse function of x^U'{x) 

D{x^U'{x)) = X. (77) 

Therefore, the mean radius is given by 



Lastly, introducing the function G(x) by 



'^ = D[^,]- (78) 



G(-) = + U{D{x)l (79) 

the binding energy of the system is given by 

i?W(rn,Q) = iVG(^). (80) 

In this case again, the D and G functions are universal since they depend on the form of the U potential, but are 
independent of the system and of the excitation number (independent of N, m, Q). 
From expression (j80l) . one deduces immediately the duality relation 



i?W(m,g) = ^i?(^')(m,|Q). (81) 



This form is exactly the same as the most general one (|23p . In this expression, the same mass appears for systems 
with TV and p particles. The duality relation leads to link between different excited states of both systems. 

Alternatively, one can use the property (|75p to obtain many other possibilities. For example, choosing the value 
(5 = N/p, one has the alternative duality relation 

E(^) (m, Q) = -i? W ( ^m, q] . (82) 



p \ p'^ y 

In this expression we decide to maintain a one to one correspondence in the spectrum but for systems with different 
particle masses. 

In ([5T]) and ([5^ . the values of N and p were fixed numbers; one can consider also p as given but N — j3p as variable. 
Thus /? is a rational variable such that I3p is an integer. Then, the last relation (|82p can be recast under an alternative 
form 

(m, Q) = I5E^P'> {(3^m, Q) . (83) 



C. Two-body interaction only 

In absence of one-body interaction, [/ = 0, we have also interesting conclusions. Let us choose the special value 
a = 2/{N - 1) in the general results of Sect. TsTKl One has VaN = N/y/C^ and Z{x) = V{x)/N. Introducing the 
mean radius tq ~ y, the transcendental equation (j69l) reduces to 



mCl, 



rlV'{r^). (84) 



Defining the D{x) and G{x) function as in (1771) and ([7^ calculated with the potential V{x), one has 



( NQ 



(86) 



N 



As the previous case studied in Sect. IIVBI the D and G functions are universal. 



13 



From expression (|86p . one deduces immediately the duality relation 

which is identical to the general case ([551) . Again, one can use the property ([75]) to obtain many other possibilities. 
For example, choosing the value /3 = y^p{N — l)/{N{p — 1)), one has the alternative simpler duality relation 

£;(^'KQ) = ^i?(-)(^™,^g). (88) 

In this last relation, let us choose /3 — -yA/V/p, one arrives at the relation 

^("),„.o,.|^E..(„.z^y|o). m 

In this expression, we decide to keep the same mass for systems with N and p particles. The duality relation leads to 
a link between different excited states of both systems. Choosing the value (3 = Cn/Cp in equation (155)) . one obtains 
the alternative expression 

In this expression we decide to maintain a one to one correspondence in the spectrum but for systems with different 
particle masses. 



D. Link between one and two-body interaction 



We consider the case of a A^-body system with particles of mass m and interacting with a two-body interaction 
V{x) only. The corresponding binding energy for a state of excitation number Q, E^^''{m, Q), is given by equation 
(|86| . Now, we consider the case of a A^-body system with particles of mass m! and interacting with one-body 
interaction of the form U{x) = cV{x) only (c being a constant). In this case, the corresponding energy for a state 
of excitation number Q', E^^\m\ Q') is given by the equation 

E^^^){m'M) = cNGi^^] (91) 



with the same definition of the function G{x). 

If one chooses the following link between the parameters. 



the arguments of the G function are the same and £'(^) is expressed as CnE'-^^ /{cN). Therefore we arrive at the 
following interesting general duality relation (using again (j75p ) 



(93) 



where /3 is any real parameter. Choosing /3 ~ 4:C\/Cn/{N — 1) , one gets 



which is a special case of ([55)) applied to a nonrelativistic system. Choosing /? = 1, one has instead 



(95) 
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giving a one-to-one correspondence between the spectra of two systems with different particle masses. Choosing 
/3 = 2y/cN /{N - 1), one has rather 

i?W(m,Q) = (m,^-^Q^ , (96) 

giving a relation between two systems with the same particle mass but for different states of their spectrum. 

V. PASSING FROM NONRELATIVISTIC TO ULTRARELATIVISTIC LIMITS 

A. General considerations 

For systems submitted to only one type of potential (either one-body potential or two-body potential) we showed 
that both the ultrarelativistic limit and the nonrelativistic limit for the eigenmasses share the property of being 
expressed in terms of universal functions. The F function for the ultrarelativistic case and the G function for the 
nonrelativistic case are independent of the system and of the excitation quantum numbers but depends only on the 
form of the potential under consideration. 

If the same form of potential is used in both situations, the expressions of the F function (see (j53p l and of the G 
function (see ([75)) ) are not identical so that the corresponding spectra are quite different. However, we will show below 
that if the potentials are different but linked by a certain relationship, one can arrive at very interesting conclusions. 

Let us start from a nonrelativistic problem which is based on a potential U{r) (no matter one or two-body). One 
has to calculate G{R) where R — R{N, m, Q) is a special combination of the parameters N, m and Q. This implies 
to calculate the x value from 

x^U'{x) = R (97) 

and then 



G{R) = ^^+U{x). (98) 



Let us consider now a new nonrelativistic problem with the same value of R but for a new potential W{r) defined by 

W{r) = C/(aVf), (99) 

where a is a dimensioned constant (in order that, in both expressions, r has a dimension of a length while U and W 
have dimension of energy), which is undetermined for the moment. 

Introducing, instead of x, the variable y defined by a: = ot^. Starting from (|97p . the y quantity results from the 
transcendental equation 

y^W\y) = S^^, (100) 



G{R) = - + W{y) = F{S). (101) 
V 



while equation (1981) is transformed into 



In this relation, the G function is the universal function based on potential U{r) whereas the F function is the universal 
function based on potential W{r). This relation allows to propose new duality relations between an ultrarelativistic 
and a nonrelativistic treatment. Let us distinguish below the case of one-body and two-body potentials. 

B. One-body interaction only 

Let us consider a system described by a nonrelativistic treatment based on a one-body potential U (r) whose binding 
energy is E^^''{m,Q). We showed in Sect. lIVBl (see equation ([50]) ) 

n2 

E^^\m,Q) ^ NG{R) with . (102) 
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If one uses rather the potential W{r) as in (|99l) . one has, due to property (|101l) . 



E'^^\m,Q) ^ NF(S) with S = —r. (103) 



If we take the special value S = Q/N, the value NF{S) represents (see ([5^ ) the ultrarelativistic mass Mi''\Q) of 
the same system but obtained with the potential W{r). The condition on S and the link between S and R leads to 
the condition defining the value of a, namely 

Q = 2mNa^. (104) 
The previous conclusions can be stated as a theorem: 

If E^^\m,Q) is the binding energy of a nonrelativistic system governed by the one-body potential U{r) 
and if Mi^\Q) is the mass of the related ultrarelativistic system governed by the one-body potential 
W{r) defined by W{r) = U { Qr / {2mN)] , then one has the general property Mi^\Q) = E'-^^m, Q). 



From its definition, the potential Wir) depends on m so that the notation M^^^ (Q) which we have used up to now, 
depends indirectly on m, as imposed by the theorem. 



C. Two-body interaction only 



Let us now consider a system described by a nonrelativistic treatment based on a two-body potential V{r) whose 
binding energy is E^^\'m, Q). We showed in Sect. IIVCI (see equation 



E^^\m,Q) = CNG{R) with i? = (105) 



If one uses rather the potential W{r) as in ([99]), one has, due to property (jlOip . 



E'^^\m,Q) ^CnF{S) with 5=-^. (106) 



If we take the special value S = 2Q/{{N — 1)^/Cn), the value CnF{S) represents (see ((6T|) ) the ultrarelativistic 
mass Mi^\Q) of the same system but obtained with the potential W{r). The condition on S and the link between 
S and R leads to the condition defining the value of a, namely 

Q = 2m^fC^c?. (107) 
In this case again, one can state the following theorem: 

If £;(^)(m, Q) is the binding energy of a nonrelativistic system governed by the two-body potential V(r) 
and \iMf \Q) is the mass of the related ultrarelativistic system governed by the two-body potential W(r) 
defined by W{r) ^ V [^J Qr / (2m^ C n)) , the n one has the general property m'^\q) ~ E^^\m, Q). 

We can mention the same remark as before concerning the indirect m dependence of Mi^\Q). 



VI. APPLICATIONS 



All the duality relations presented in the previous sections are exact for the AFM solutions of quantum systems. 
We can wonder to what extent these constraints are satisfied for the corresponding exact solutions. In this section, 
we will examine the relevance of some duality relations for genuine solutions of particular systems. 

We will not test all the duality relations which are presented in this paper (some of them are consequences of others 
and do not need a second check), but we will look at the cases of two-body interactions only {U = 0) and one-body 
interactions only {V = 0), as well as nonrelativistic and ultrarelativistic kinematics. In a first study, we will consider 
nonrelativistic systems, for which relation ([75)) is specially important and will be used intensively. 
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A. Nonrelativistic systems 

1. General considerations 

Let US call e'^^(m; {rii, Zi}) the exact eigenenergy of the {rii^li} = {rii, /i, n2, Z2, ■ • ■ , ^T-Ar-i, ^w-i} state of the 
Hamiltonian 

N 2 ^ 

^ = E|^+ E nK:-r,|), (108) 

1=1 ■i<j=i 

while e{m;n,l) = e^'^\m;{n,l}) gives the exact spectrum of the corresponding two-body problem. E^^^{m;Q^-^^) 
and E{m; Q) are the approximate AFM eigenenergies of the same states. The duality relations are exactly true 
for the E values but not for the e values in general. In this Sect. IVI A[ we will consider only systems of boson- like 
particles, in order that the ground state is given by {ni,li} = {0,0}. 

Let us assume that we are able to obtain, for each set of parameters {m,n,l), the value e{m;n,l) for a two-body 
system. Nowadays, it is not difficult to get numerical solutions for such a problem. We can compute the corresponding 
value E{m; Q) analytically or, if not, numerically. Equating both values defines the Q value that leads to the exact 
value in the AFM expression. Of course, this value depends on (n, I) but also, most of the time, on m. We note this 
value Q(rn] n,l) which is thus defined by the formal equality 

E{m;Q{m;n,l)) = e{m;n,l). (109) 

For the harmonic oscillator (HO), the Q value is m-independent and is equal to 2n + Z + 3/2, whereas for the Coulomb 
potential the Q value is still m-independent and is equal to n + I + I. We showed that for a lot of potentials, the 
rn-dependence of Q is not crucial but that a good dependence in (n, I) is of major importance [T^. For the rest of this 
section, we suppose that we adopt a form of Q depending on n and Z, but not on m, which gives quite good results. 
Consequently the relation (|109p does not remain an equality but holds only approximately 

E{m- Q{n, I)) « e(m; n, I). (110) 

The duality relation ([75]) allows to express the AFM energy in terms of the energy of the ground state but for a 
different mass. Denoting (5(0,0) = Q2, one has explicitly 

E{m;Q{n,l)) ^ E{m{m;n,l);Q2), (111) 

with the definition of the fh(rn; n,l) mass 

These equations are equalities. Using them, with the approximate relation (jllOl) . one arrives at the approximate 
duality condition concerning the exact states 

e(m; n, I) « e{fh{'m; n, I)] 0, 0). (113) 

Let us define the / function by 

/(m) = e(m;0,0). (114) 

This function is universal in the sense that it depends only on the form of the potential V{r) and it can be computed 
once for all. The relation pi3p can then be recast under the form 

e{m; n, I) w f{rh{m; n, I)). (115) 

The conclusion is very strong. It means that the entire spectrum of a given two-body system can be obtained, at 
least approximatively, by using in an universal function, /(m), the various arguments m{m;n,l) given by (|112p . The 
quality of the results depends of course on the ability for the function Q{n, I) to reproduce the exact results in a 
satisfactory way. The test of duality relation (jllSp in a realistic case is the subject of Sect. IVI A 21 

The relation (|75|) being valid for an arbitrary number of particles, all the demonstrations that we have developed 
above for the two-body problem apply as well for the iV-body problem. One can define formally a principal quantum 
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number Q^^\m; {n^, li}) through an equahty similar to (I109p . The crucial hypothesis is that this value depends only 
slightly on the m parameter and can be very well approximated by a simpler function Q'''^^ ({r^i, Denoting this 
time the principal quantum number for the ground state as Q^v = Q^^\{0, 0}) the equivalent of equation (lllip now 
writes 

i;W(m;QW({n„Z,;})) =£;(m(m;{n„/J);QAr), (116) 
with the definition of the m{m] {ui, li}) mass 

-K{".,^.})=(Q(^)(^^^^y)),m. (117) 

Thus, one has an approximate duality relation concerning the A^-body spectrum in term of the ground state of the 
same system for another mass 

eW(m; {n„ k}) « eW(m(m; {n„ k}); {0, 0}). (118) 

The test of duality relation (|118p in a realistic case for iV = 3 will be discussed in Sect. IVI A 31 

Finally, a link can be found between the A^-body and the 2-body systems. Since a link between excited states and 
the ground state has been proposed above, it is sufficient to search for a relationship between ground states. In order 
to do that, let us apply relation ([55)1 for p ^ 2 leading to (with again the use of ([75]) ) 

E'^^\m; Q) = CNE{2m/N; Q/Cn) = C n £{20^ m/N; uQ/Cn), (119) 

where a is an arbitrary real parameter. Let us choose the value m = ffi{m] {rii, ?.;}), Q = and a ~ CnQ2/Qn in 
the previous equation. One gets 

(rfi{m; {n„ /J); Qn) = CNE{M{m; {n,, k}); Q2), (120) 

with the definition of the M mass 

M(m; {n., / J) = ^ ( ' . (121) 



If the AFM energies give a good approximation of the exact results, one can expect that i?(^)(m(m; {rii, /i}); Qat) « 
eW(m(TO;{n„?,});{0,0}) and E{M{m; {n„k}); Q2) « e(M(TO; {n„ ?,}); 0, 0) = f{M{m;{n,,k}). Owing to the 
relations (jllSp and (I120p . one arrives at the very important result 



e 



W(m; {n„ k}) « Cjv/(M(m; {n„ k}). (122) 



The conclusion of this relation is even stronger than (jllSp . Equation (|122p proves that the whole spectrum of all 
systems can be obtained approximatively by the calculation of a universal function, /(m), corresponding to the 
ground state of the 2-body system with the same potential, and for arguments M(m; {rii, k}) given by (|12ip . Getting 
/(m) is a very easy task. Solving the 2-body system can be performed with a great accuracy for any potential (let 
us recall that this potential must not depend on m and N)] moreover, obtaining the ground state energy is free from 
possible numerical complications arising for excited states. Let us note also that the duality relation (jl22p . which 
is obviously an approximation, concern the exact eigenvalues only; the AFM values, which were very convenient 
intermediate quantities in our demonstration, have completely disappeared. Testing the link between the ground 
states of a 2-body and a 3-body system is the subject of Sect. IVI A 41 

Let us point out an additional comment. All what we did for the ground state of a two-body system can be 
reproduced identically for any other excited state (no,/o)- It is possible to introduce a universal function /no,io('^) 
built on this state. Then, (jl22p can be expressed in term of the fno,io function, but with a new argument given by 
(|121l) in which Q2 has been replaced by Q{no, lo)- 

The scaling laws or the dimensional analysis allow generally to rewrite a Hamiltonian on a reduced form, easier 
to study. The eigenenergies can then be expressed in terms of a dimensioned parameter (an energy scale) and other 
dimensionless parameters which depends on the physical quantities of the problem (particles masses, length scales, 
etc.) but not on the quantum numbers of the corresponding eigenstates. That is why equations (jllSp and (jl22p are 
both of different nature than the usual scaling laws. 
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2. Testing a two-body system 

In this section, we test the duahty relation (jll5p for two particles interacting with a pure linear potential V{r) = r, 
so that the corresponding Hamiltonian in reduced variables is simply 

H=— + r. (123) 

TO 

The first thing to do is to obtain the ground state energy as a function of the mass in order to plot the universal 
function /(to) for the linear potential. The numerical value has been calculated with a great accuracy using the 
so-called Lagrange mesh method [l^. The form of the /(to) function is shown in Figure [TJ This function decreases 
very rapidly from infinity at the limit to — >■ to around 2 for to = 1; then, it decreases very slowly to zero for large 
values of to. 

12 

10 - 

8 - 

/( ) 





1 2 3 4 5 




FIG. 1: Universal function /(m) as a function of the mass m for Hamiltonian (|123|l . 



Now we calculate, still using the Lagrange mesh method, the spectrum of a system with mass m = 4, in order to 
match the results of @. We think that limiting ourselves to states with < n < 3 and < / < 3 is enough for our 
purpose; the corresponding eigenenergies e{m = 4; n, I) are given in the first line of Table HI 

In order to calculate the effective masses rh{A;n,l) as defined in (jll2p . we need a prescription for the choice of 
the principal quantum number. We study here two possibilities: the HO prescription Q{n,l) — 2n + Z + 3/2 and 
the improved value proposed in [5j Q{n,l) = 1.789n + I + 1.375. For both prescriptions, we calculate the effective 
masses (jll2p and insert them in the /(to) function to obtain values f{rh{A]n,l)) which are ultimately compared to 
the exact results e(4;n, /). The approximate values f{'fh{A;n,l)) are reported in the second line of Table U for the 
HO prescription whereas the third line of the same Table shows the results with the improved formula. Note the the 
ground state is always exactly computed by definition. 



TABLE I: Values of e(m = 4;n,l) for various approximations. For each set {n,l), the first line is the result obtained by 
numerical integration e(4;n, Z), the second and third lines are the results coming from duality relation /(m(4;n, Z)) with the 
HO prescription Q{n,l) = 2n + I + 3/2 and the improved prescription Q{n,l) — 1.789n + / + 1.375, respectively. Between 
parenthesis is also indicated the deviation (%) of the approximate result as compared to the exact ones. 

I e{4;0,l) e(4;l,0 e(4; 2, Q £(4; 3, /) ~ 

1.473 2.575 3.478 4.275 
1.473 (0) 2.591 (0.6) 3.502 (0.7) 4.307 (0.7) 
1.473 (0) 2.567 (0.3) 3.461 (0.5) 4.251 (0.5) 

1 2.117 3.077 3.911 4.665 
2.071 (2.2) 3.064 (0.4) 3.915 (0.1) 4.682 (0.3) 
2.120 (0.1) 3.083 (0.2) 3.913 (0.05) 4.662 (0.06) 

2 2.676 3.546 4.327 5.046 
2.591 (3.2) 3.502 (1.2) 4.307 (0.5) 5.042 (0.08) 
2.680 (1.5) 3.559 (0.4) 4.339 (0.3) 5.055 (0.2) 

3 3.182 3.989 4.728 5.416 
3.064 (3.7) 3.915 (1.8) 4.682 (1.0) 5.390 (0.5) 
3.186 (0.1) 4.005 (0.4) 4.746 (0.4) 5.434 (0.3) 



From the table one sees that the duality relation (jllSp is very well satisfied in any case. The very simple HO 
prescription is less good (specially for large I values) but remains acceptable with discrepancies around 1%. The 
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improved prescription, which raises the degeneracy, is very impressive (even for large I values) with an average 
precision of the order of 0.5%. 

We checked also the validity of this duality relation for other types of potential (square root [vr^ + a], funnel 
[~a/r + br]) and find that the accuracy is always of the order of 1% or better. Thus, we think that the studied duality 
relation is well adapted for the spectrum of a two-body problem. 



3. Testing a three-body system 

In this section, we test the duality relation (IllSp for three particles, with a mass m, interacting with a pure linear 
potential V{r) = r, so that the corresponding Hamiltonian in reduced variables is simply 

^ = Eil+ E l-^-^.l- (124) 

i—1 i<j — l 



To compute the solutions of this Hamiltonian, we use a variational method relying on the expansion of trial states 
with a HO basis [16,]. We can write 

S=0 q{B) 

where B = 2(ni + 712) + (^i + I2) characterizes the band number of the basis state and where q summarizes all the 
quantum numbers of the state (which can depend on B). This procedure is s pec iall y interesting since the eigenstates 
of H are expanded in terms of AFM eigenstates (up to a length scale factor) [111, [ij]- In practice, a relative accuracy 
better than 10^^ is reached with Bmax = 20. Such results are denoted "exact" in the following. Some exact eigenvalues 
e of (jl24p for m = 2 are presented in Table HIl with the quantum numbers {ni, Zi, n2, h} and the B value of the main 
component in the expansion (jl25l) . 

In order to compute the effective mass (I117|) . we need a prescription for the principal quantum number. We 
consider two formulas. The first one is the HO result: Qho = 2(ni + 712) + (/i + 12) + 3. The second one is given by 
QwKB — "^l*^! + "2) + (^1 + h) + 3. This prescription comes from the fact that the ratio between the coefficients 



of n and / for a nonrelativistic particle in a linear potential is predicted to be by a WKB method 17 [ and the 

AFM 5]. The set of quantum numbers used for a state is fixed by the quantum numbers of the main component 
in the expansion (jl25p . It is then possible to compute the energy of an excited state from (jll7p and (IllSp . The 
corresponding values denoted e((3Ho) and e(QwKB) are presented in table HIl Note the the ground state is always 
exactly computed by definition. 



TABLE IL Some exact eigenvalues e of p24p for m — 2 with the quantum numbers {ni, Zi, n2, ^2} and the B value of the main 
component in expansion (|125p . The brackets indicate that two components are equally present with ni -f^ n2 and h ^ l2. 
Energies predicted by (|117p and (|118p are given in columns e(QHo) and e(QwKs) following the prescription chosen for the 
principal quantum numbers (see text). Between parenthesis is also indicated the deviation (%) of the approximate result as 
compared to the exact ones. 
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(0 


5) 


3 


[1,1,0,0] 


7.608 


7.726 


[1 


6) 


7.566 


(0 


6) 




1,0,0,1 


7.702 


7.726 


'fi 


3) 


7.566 


(1 


8) 




[0,2,0,1] 


7.854 


7.726 


[1 


6) 


7.726 


(1 


6) 


4 


[2,0,0,0] 


8.309 


8.562 


J3 


0) 


8.256 


(0 


6) 




1,1,0,1 


8.391 


8.562 


[2 


0) 


8.410 


(0 


2) 




[1,2,0,0] 


8.426 


8.562 


[I 


6) 


8.410 


(0 


2) 




0,1,1,1 


8.572 


8.562 


'fi 


1) 


8.410 


(1 


9) 




0,2,0,2 


8.707 


8.562 


[1 


7) 


8.562 


(1 


7) 
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From this table one sees that the duahty relation (|117l) - (|118p is very well satisfied in any case. The very simple 
Qho prescription already gives good results, but the Qwkb prescription raises partly the degeneracy and improves 
appreciably the global agreement. 



4. From two- to three-body ground states 



Let us consider a system of particles interacting via two-body forces and let us assume that the principal quantum 
number of the ground state can be written Qn = p{N — 1) with p independent of N. It is true for the A^-body HO with 
p = 3/2 (see ([6|)). In this case, using (jl20l) and (I121|) . and assuming that the AFM solution is a good approximation 
of the true solution, we can write 

e'^)(m;{0,0}) « CNe{Nm/2;0,0). (126) 

We have then an approximate link between the ground state of N particles interacting via a two-body potential and 
the ground state of the corresponding two-body system. We have checked that, for = 3, relation (jl26p is satisfied 
with an error less that 1% for V{x) = ax, around 1% for V{x) = —a/x -\- bx and around 6% for V{x) = —a/x (these 
potentials deviate more and more from a quadratic one). 

The three tests performed in sections IVI A21 IVI A 31 and IVI A 41 prove that the most general formulation of the 
duality relations, as expressed by (jl22p . is indeed satisfied with a rather good accuracy of the order of a few percent, 
at least for = 2, 3. 



B. An ultrarelativistic system 



Let us consider the following ultrarelativistic Hamiltonian 

N N 

H = Y.^,+aY,\n-Rl (127) 

1=1 i=l 

which presents some interest for hadronic physics [l3| and whose AFM eigenmasses are given 

In [ll| , it has been shown that the accuracy of this formula can be improved as compared to the pure HO prescription 
in the case iV = 3 by the determination in ^ of: 

• The 7 term by using an accurate variational solution for the ground state. 

• The ratio ai/Pi by using the characteristics of a WKB solution for an ultrarelativistic particle in a linear 
potential 17]. 

The same method can be used in the case N = 2. Finally, we obtain 

M^2) ^ ^8aQ(2) with q(^) ^ln + l + -, (128) 

2 TT 

Mi""^ ^ J—aQ(^) with Q^^^ ^^{ni+n2) + h+l2 + 3. (129) 

V TT Z 

The relative errors are below 2% for (I128p and below 1% for (|129p (at least for quantum numbers such that n < 3 
and / < 3). So, we can consider that these formulas are very close to the exact solutions. 

We can now check if relation (|55p is, at least approximately, verified with iV = 3 and p = 2 for the (nearly) exact 
solutions of this system. A priori, we could expect that Ali'^\Q^^'>) « 3/2 M^^^ (2(3(^^/3). The calculation gives 



This must be compared with (I129D . Since 32/7r w 10.2, the relative error is around 15%. 
But, we can also hope that Afl^'(Q(2)) ^ 2/3 Afi^^(3Q(2)/2). The calculation gives 



This must be compared with (jl28p . Since 64/ (37r) « 6.8, the relative error is also around 15%. In these cases, we can 
consider that the relevance of the duality relation is demonstrated. 
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C. An ultrarelativistic-nonrelativistic duality 



The eigensolutions of these two two-body Hamihonians 



m 



br. 



are given respectively by 



EniQn) — 




„, with Qn = 2n + I + -, 
m 2 



MuiQu) ~ V^bQu with Q„ = -n + / + -. 

2 TT 



(132) 
(133) 

(134) 
(135) 



Relation (jl34l) is exact while the relative accuracy of (|135|) is around 1-2%. If we impose Qn — Qu — Q, the duality 
relation presented in Sect. |V]is verified by setting 2bm — aQ as shown above. 

Let us examine what becomes this property for the genuine solutions. If we assume that En is known, we can 
compute the exact solutions of with b = aQn/{2m). For n < 3 and I < 3, the relative error compared to (|135p 
is below 10%. On the contrary, if we assume that M„ is known, we can compute the exact solutions of Hn with 
a = 2bm/Qu- For n < 3 and I < 3, the relative error compared to (|134p can then reach 38%. Following the physical 
situation considered, the accuracy of this kind of duality relation can strongly vary. But, in some cases, the error 
could be quite small. 



VII. CONCLUSIONS 



In this paper, we investigate the possible relationships between the energies of particular states for a given system 
and the energies of other states for another system. We call these relationships "duality relations" . Systems with 
arbitrary number of particles are considered but in a well defined framework. The basic hypotheses are the following 
ones: 

• All the particles of the studied systems are identical. 

• We do not take into account internal degrees of freedom for particles and, for convenience, we focus only on 
bosons in applications (the case of fermions can be dealt with as well, but it is more involved and is not studied 
in this work). 

• We consider the kinetic energy operator either of a nonrelativistic form or of a semirelativistic form. 

• We limit ourselves to one-body and two-body types of potentials. 

• The potentials do not lead to coupled channel equations, so that we are faced with a Schrodinger equation for 
a nonrelativistic kinetic energy, and with a spinless Salpeter equation for a semirelativistic kinetic energy. 

• The potentials entering the formalism do not depend on the mass m nor on the number of particles N. 

Despite the limitations, the number of systems which can be concerned by this work is quite large and may represent 
a lot of different physical situations. 

As an intermediate tool, we introduce the auxiliary field method (AFM) to get approximate values of the eigenener- 
gies. The important parameters are the mass m and the principal quantum number Q{{ni, U}), function of the various 
radial quantum numbers Ui and orbital quantum numbers li. The form of this function is not given by the theory 
and is left to the cleverness of the physicist. Usually, the harmonic oscillator prescription ([5]) is a good starting point, 
but a more elaborated formula, such as the one given by ([7]), can improve substantially the quality of the results. 
The AFM provides an expression £;W(m;Q) of the eigenenergies (or eigenmasses for a semirelativistic kinematics) 
in terms of these parameters. In a series of papers |5l4l2l|. we proved that, the AFM approximation can be calculated 
analytically for many systems and leads to encouraging results. 

A duality relation is, at the beginning, a link between AFM approximations E'^^\m; Q) and i?(P^(m'; Q'). The case 
p = 2 is particularly important because the corresponding duality relation allows to obtain the spectrum of a iV-body 
system in terms of the spectrum of a 2-body system, a much more favorable situation. We obtained very general 
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duality relations when only one-body or two-body potentials are present. More interesting conclusions can be drawn 
in the ultrarelativistic limit (m = 0) or in the nonrelativistic limit. In these particular cases, we showed that the AFM 
results can be expressed in term of a unique universal function for a given potential and for a given kinematics. This 
crucial feature has a lot of important consequences for the duality relations. We proved for example that the spectrum 
of a nonrelativistic system governed by a two-body potential (the same property holds for a one-body potential) V{r) 
is the same as the spectrum of an ultrarelativistic system governed by a two-body potential W{r) = V{a\/r) provided 
that the parameter a is chosen adequately. The nonrelativistic kinematics allows additional sympathetic features. In 
particular, it was shown that the spectrum of such a iV-body system can be obtained from the ground state of the 
corresponding two-body system. 

The duality relations presented here are exact for the AFM approximation of the eigenenergies. Assuming that this 
approximation leads to results close to the exact ones, the duality relations can also be applied to the exact states. 
This means that we have an approximate link between the exact energies of a A^-body system e^-^\m; {n^, k}) and 
the exact energies of a p-body system e'^P^m'; {n-, /■}). The case p = 2 is specially interesting because nowadays we 
are able to get easily the spectrum of a two-body system numerically with a great accuracy. 

In particular, for a nonrelativistic system, we obtained a very important conclusion. The whole spectrum of any 
system (arbitrary N value) can be calculated, approximatively, through a universal function /(m) which is nothing 
else than the dependence of the ground state of the two-body system (with the same potential) on the mass. This 
function is very easy to obtain numerically in any circumstances. The duality relation is expressed by equation (jl22l) . 
It is sufficient for obtaining the searched value e^^^{m;{ni,li}) to calculate the universal function for an argument 
M(m; {n„li}) given by 

The important point that we want to stress is the following one. The duality relations are exact for the AFM 
approximations of the eigenenergies while they are only approximately true for the exact eigenenergies (the quality 
probably deteriorates with increasing values of N). This drawback is compensated by the fact that the corresponding 
duality relations concern the exact eigenenergies without any reference to the AFM expressions which were only 
an intermediate tool. Thus, the consequences of duality relations would remain true independently of the AFM 
approximations. We propose to apply all the duality relations presented above for the exact eigenvalues, but only in 
an approximate way. 

The duality relations have been tested in a number of interesting cases. In particular, using a linear potential for 
N = 2 and 3, we showed that they are fulfilled with an accuracy of the order of 1-4%. For other types of potentials, 
for larger values of N and quantum numbers {ni,li}, the quality is likely poorer, but we think that these duality 
relations are nevertheless able to give valuable informations concerning complicated systems which are not easy to be 
dealt with, analytically or numerically. 



Appendix A: Tests of duality relations for AFM energies 



In order to test the duality relations presented above, we can use the semirelativistic A'-body harmonic oscillator 
studied in [III 



N 



N 



R- 



N 

E 



(Al) 



2—1 i=l i-cij — l 

The AFM solutions, for the general, the ultrarelativistic and the nonrelativistic cases respectively, are given by 



2Nm 



VSYn 
3 



1 



G-{Yn) 



+ G_(rjv)' 



with Yn = 



2N' 



{k + pN)Q-- 



2/3 



^ ±[2N{k + pN)Q^]"\ 



-{k + pN)Q. 



(A2) 
(A3) 
(A4) 



Let us mention that the last formula is an exact solution. The function G_ (Y) is the root of a quartic polynomial 
whose explicit form is given in [10] for instance. Using these formulas, it is easy to check the relations (P^. psp . 

(ESI), (IMI), (Uni), ([23, (EHi-dlsl), (EZli-dMl), dMli-dMl). 

The semirelativistic harmonic oscillator with one variable, studied in [10] , has the following form 



(A5) 
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The AFM solutions, for the general and the ultrarelativistic cases respectively, are given by 



2(TTO 



f 16a- \ 



(A6) 
(A7) 



Using these formulas, it is easy to check the supplementary relations ((57)l . (1551) . ([M)) . 

To check the duality relations between ultrarelativistic and nonrelativistic systems with an example, we can consider 
the following Hamiltonian 



N 



N N 

sgn(A)a^|r, -Jl|^+sgn(A)6 ^ In-r^l 

1=1 i<j = i 



(A8) 



Defining the quantities 



^A = a|A| ( - 



and Bx = b\X\N 



N -1 
2Q 



the AFM solutions, for the nonrelativistic and the ultrarelativistic cases respectively, are given by [11| 



A + 2 

2A 
A + 1 



Q 



[Ax + Bx 



A 



(A9) 

(AlO) 
(All) 



Using these relations with potentials W{x) cx and U{x) cx x^^, the two theorems presented in Sect. |V] can be 
checked. 
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